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Abstract 

Seabirds are affected by changes in the marine ecosystem. The influence of climatic factors on marine food webs can be 
reflected in long-term seabird population changes. We modelled the survival and recruitment of the Mediterranean storm 
petrel {Hydrobates pelagicus melitensis) using a 21 -year mark-recapture dataset involving almost 5000 birds. We 
demonstrated a strong influence of prebreeding climatic conditions on recruitment age and of rainfall and breeding period 
conditions on juvenile survival. The results suggest that the juvenile survival rate of the Mediterranean subspecies may not 
be negatively affected by the predicted features of climate change, i.e., warmer summers and lower rainfall. Based on 
considerations of winter conditions in different parts of the Mediterranean, we were able to draw inferences about the 
wintering areas of the species for the first time. 
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Introduction 

In the face of the overwhelming data on the effects of climate 
change, it is crucial to understand how wildlife will react to such 
changes. Seabirds are considered sentinels of environmental 
change [1], as they range over large areas of the ocean and are 
therefore directiy affected by any changes that may occur 
[1,2,3,4]. Climate change affects oceanographic conditions and 
the availability of seabirds' prey, ultimately influencing population 
dynamics [5,6,7,8]. Changes in fluid circulation may also induce 
positive consequences as it has been evidenced in the foraging 
efficiency and breeding success of albatrosses due to changes in 
wind regimes [9]. 

In view of the effects of environmental variability on the prey of 
seabirds, the climate-induced warming of oceanic surface waters to 
levels that exceed the boundaries of a fish species' thermal niche 
can have a direct effect on the distribution of the species. Even 
within this thermal niche, subtle (<0.5°C) temperature changes 
can have profound effects on growth, survival and reproduction 
[1,10]. As a consequence, climate change can have an indirect 
impact on seabirds [1 1] influencing prey quality in marine food 
webs [12,13] and complex bottom-up and top-down processes in 
ecosystems [14]. Climate factors may affect breeders directly 
through changing the distribution of important prey species or 
indirectly through affecting spawning and recruitment of impor- 
tant fish prey species giving lagged effects on seabird's demo- 
graphic traits. Changes in prey abundance and quality during pre- 
breeding and chick rearing period may result in variation in adult 
body condition affecting their quality as breeders and or their 
survival. In that case changes that result in reduced caloric value of 
prey, compared to those occurring under "normal" conditions, 



may affect reproduction and chick survival [15,16]. Chick survival 
may be indirectly affected by the same factors, influencing growth 
rate and body condition at fledging, but also by colony conditions 
in terms of humidity influencing thermo-regulation. 

Long chick-rearing periods involving both parents are a 
characteristic of seabirds [17]. Reproductive output may therefore 
be strongly impacted by adverse conditions that force the birds to 
fast for longer periods or to make excessive metabolic effort [18], 
as a consequence juvenile survival and the overall population 
trend are at risk. 

Cyclic and extreme climatic events are inherent to oceanic 
systems and have occurred long before human-induced global 
warming. Seabirds have evolved in this fluctuating environment 
and are negatively affected by the rapid increase of environmental 
stochasticity as a consequence of global warming [19]. 

In the Mediterranean, models have predicted an increase in 
temperature and a decrease in rainfall [20], which could produce 
the effects mentioned above, affecting the whole food web. In this 
context, we aimed at linking the effect of climatic conditions to the 
demography of seabirds with particular reference to the Mediter- 
ranean storm petrel (Hydrobates pelagicus melitensis). Our aims were: 

1 . to determine age of recruitment in Mediterranean storm petrel, 

2. to investigate effects of environmental conditions in different 
time of the year on important demographic traits and 3. to reveal 
important wintering-area for this population. The last point is 
particularly challenging in species whose distribution is unknown 
and where technology is not yet available to track small animals at 
sea beyond the breeding period. 

The Mediterranean storm petrel is a taxon distinguished both 
morphologically and genetically from the Adantic populations 
[21,22]. On the IUCN Red List the storm petrel is not considered 
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to be under any threat and is thus categorized as Least Concern 
[23], but if we focus only on the Mediterranean subspecies, it may 
be endangered mainly due to habitat degradation and introduc- 
tion of predators, such as rats and cats, requiring conservation 
awareness [24]. Breeding distribution and abundance is probably 
under estimated, due to the highly pelagic attitude and nocturnal 
behaviour [25], The known Mediterranean population is estimat- 
ed at between 11 000-16 000 breeding pairs, with three main 
population cores identified in Malta, Sicily and the Balearic 
Islands [26] . The current trend of the Mediterranean storm petrel 
populations is stable [26], although there is a lack of basin wide 
studies. In the few sites where mid-long term monitoring has taken 
place (Marettimo, Filfla and Benidorm), it has been observed that 
the colonies remain stable. 

Dietary needs have evolved according to Mediterranean basin 
food-web and are mainly supplied by small pelagic fish [27]. 
Considering fluctuations of small pelagic fish whose biomass 
estimates resulted to be negatively correlated with the mean sea 
surface temperature [28,29,30,31] we predict that climatic 
variation both in sea surface temperature and in environmental 
temperature and rainfall may affect storm petrel survival at 
different stages of their life cycle. Our hypothesis is that conditions 
experienced at fledging can influence first year survival and 
breeders may be affected by environmental conditions during the 
pre-breeding and breeding period with effects on age of 
recruitment and/or adult survival. We used temperature and 
rainfall in the pre-breeding and breeding period and sea surface 
temperature during winter in testing whether and in what measure 
birds are affected directly, chlorophyll concentration was used as 
indirect covariate. 

We investigated the effects of climatic conditions on the 
demography of this taxon by analysing the life histories of almost 
5000 individuals monitored since 1991 in the Egadi Archipelago, 
Italy. We then inferred that local conditions from the different 
areas during the winter can affect vital rates and analysed 
correlations between trends in different Mediterranean upwelling 
zones and storm petrel demography in order to predict the 
wintering areas. 

Materials and Methods 

Study area and species 

The storm petrel {Hjdrobates pelagkus) is the smallest procellarii- 
form present in the Mediterranean Sea, where the subspecies 
melitensis lives. This species lays one egg in May-June and chicks 
fledge in late August-September. Mark-recapture data were 
obtained from Marettimo Is. (37°58'N, 12°03'E, Figure SI), Italy, 
between 1991 and 2012. The main colony, placed in a large cave 
accessible only from the sea, is the only accessible colony consisting 
of approximately 2500 pairs [25]. Recruitment in Marettimo is 
observed earlier than in other colonies [32,33]. Each year, the 
colony was visited at least once between June and August, and 
adults and chicks were captured at nest and banded with metallic 
inox rings for a total of 4926 marked individuals. Sexing of the 
birds began in 2007 on chicks and adults following the methods 
already described in Albores et al. 2010 [34]. Adults caught at the 
colony are all breeders, as pre-breeders and skippers are not 
present. 

Climatic covariates 

We used two types of covariates: time-dependent individual 
covariates and environmental covariates (Table 1). Individually 
based covariates follow individuals life history with no relation 
with other individuals while values of environmental covariates are 



shared by all individuals with no reference to their life history. In 
the case of time-dependent individual covariates, we had one value 
per individual, referring to the environmental conditions experi- 
enced at birth and the subsequent covariate values will charac- 
terize individual life history: the sequence of covariate values for 
each individual starts in the year of birth, hence birds born in 
different years have different sequences for the same covariate 
(e.g.: WNAO for the first 5 years of life of a bird born in 1995 will 
be: 1,363; -0,620; -0,067; -0,227; 0,643; while for a bird born 
in 1998 WNAO values will be: -0,227; 0,643; 1,303; 0,040; 
0,237). We assumed these covariates to characterize individuals' 
life history (as for example in the case of individuals' weight), 
although they are not dependent on an individual's life history. In 
the case of environmental covariates, we had a single value per 
year that was shared by all individuals. Integrative indexes used as 
individual and external covariates were obtained performing 
Principal Components Analysis, PCA, see below for details on 
variables used [35]. 

We defined four time-dependent individual covariates (WSR, 
western Sicily rainfall, WNAO, winter NAO, PB, pre-breeding 
period integrative covariate, and B, breeding period integrative 
covariate) to test climate effects on the recruitment probability and 
survival of chicks and adults. 

• WNAO: we selected the December through February winter 
North Atlantic Oscillation [36] due to its broad geographic 
coverage and used the average value of the selected period. 

• WSR: to define the integrative covariate a PCA was performed 
on rainfall data from the four meteorological stations closest to 
the island of Marettimo: Erice, Trapani, Marsala, and Mazara 
del Vallo (Figure SI) [37,38]. The analysis considered the total 
rainfall recorded in 12 months relative to the life cycle year 
(May to April, considering the year starting with egg laying). 

• PB and B: to incorporate the uncertainty with which climatic 
variables affect the pre-breeding (March-April) and breeding 
(May-August) periods, we included temperature, rainfall, and 
NAO anomalies in the PCAs to calculate the integrative 
covariates as Pre-Breeding, PB, and Breeding, B, indexes. 

To test winter period climate effects on adult survival and 
identify potential wintering sites for the species, we considered two 
external covariates defined as environmental covariates: SST, sea 
surface temperature, and CHL, chlorophyll concentration. We 
considered time series from 2001 to 2010 of SST and CHL as 
proxies of environmental conditions [39] from five areas with 
upwelling conditions in different regions of the Mediterranean Sea 
(Alboran Sea, Balearic Sea, Sicilian Channel, Adriatic Sea and 
Aegean Sea) [40]. To test the effect of winter conditions, we 
analysed data from December to February, performing 2 PCAs for 
each sea area (1 PCA for sea surface temperature, SST, and 1 
PCA for chlorophyll concentration, CHL, each one based on 5 
points). We then used the PC factors as external covariates in the 
model for a total of 1 1 covariates tested. 

Survival and recruitment estimate from encounter history 

We applied a multistate mark-recapture analysis to estimate 
chick and adult survival, to assess recruitment and to evaluate the 
potential effects of time and sex on these variables [41,42]. 
Recruitment probability is given by the probability of transition 
from the non-breeder to the breeder state and can vary according 
to age and to time. We constrained variation with age by 
considering that recruitment probability stabilizes to a constant 
value from a given age onward. We defined first year survival after 
fledging as chicks survival and adults survival that of adults of 
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Table 1. All covariates used in the analyses, and their main characteristics. 





Kind of covariateName 


Characteristics 


label 


Reference period 


Available years 


demographic 
parameters that could 
be influenced 


Individual 


Winter North Atlantic 
oscillation 


global-scale covariates 


WNAO 


December-February 


1991 


-2012 


Survival 


Individual 


Temperature anomaly. 
Pre-Breeding period 


local scale breeding ground 
covariate 


PBT 


March-April 


1 991 


-2012 


Survival and breeding 


Individual 


Rainfall anomaly. Pre- 
Breeding period 


local scale breeding ground 
covariate 


PBR 


March-April 


1 991 


-2012 


Survival and breeding 


Individual 


North Atlantic 
oscillation. Pre-Breeding 
period 


global-scale covariates 


PBNAO 


March-April 


1 991 


-2012 


Survival and breeding 


Individual 


Integrative climatic 
index 1 Pre-Breeding 
period 


PCA derived integrative 
covariates 


PC1PB 


March-April 


1 991 


-2012 


Survival and breeding 


Individual 


Integrative climatic 
index 2 Pre-Breeding 
period 


PCA derived integrative 
covariates 


PC2PB 


March-April 


1991 


-2012 


Survival and breeding 


Individual 


Temperature anomaly. 
Breeding period 


local scale breeding ground 
covariate 


BT 


May-August 


1991 


-2012 


Survival and breeding 


Individual 


Rainfall anomaly. 
Breeding period 


local scale breeding ground 
covariate 


BR 


May-August 


1991 


-2012 


Survival and breeding 


Individual 


North Atlantic 
oscillation. Breeding 
period 


global-scale covariates 


BNAO 


May-August 


1991 


-2012 


Survival and breeding 


Individual 


Integrative climatic 
index 1 Breeding region 


PCA derived integrative 
covariates 


PC1B 


May-August 


1991 


-2012 


Survival and breeding 


Individual 


Integrative climatic 
index 2 Breeding region 


PCA derived integrative 
covariates 


PC2B 


May-August 


1991 


-2012 


Survival and breeding 


Individual 


West Sicily Rainfall index 


PCA derived integrative 
covariates 


WSR 


May-April 


1991 


-2012 


Survival and breeding 


External 


Sea surface temperature 


PCA derived integrative 
covariates 


SST 


December - February 


2001 


-2010 


Survival 


External 


Chlorophille concentration 


PCA derived integrative 
covariates 


CHL 


December - February 


2001 


-2010 


Survival 



In the table we reported time-dependent individual covariates as "individual" and environmental covariates as "external". 
doi:1 0.1 371 /journal.pone.0094526.t001 



subsequent age classes. The assumptions underlying our model 
are: a. detection probability of non-breeders and skippers is zero; 
b. once in the breeder state, individuals cannot go back to the non 
breeder state; c. there is an age from which the probability of 
becoming a breeder (probability of transition from non breeder to 
breeder) levels off at a constant value. As a consequence, at a given 
sampling occasion, an animal may be a pre-breeder (a juvenile of 
one year or more, not yet recruiting, state PB), breeder (adults after 
first breeding are generally considered breeders, state B), or dead 
(state dead, unobservable). The following observations may be 
made: '1' (if detected as chick), '2' (if detected as adult) and '0' (if 
not detected). We tested the full recruitment age (the age at which 
any bird has become a breeder) and the age of full breeding 
(considering the probability of becoming a breeder as a function of 
age). To address age variation in the probability of becoming a 
breeder, we used only individuals of known age, i.e., only those 
ringed as chicks (3135), considering ringing as chicks as first 
capture event. We define the initial state vector II, the transition 
matrix cD v F (survival multiplied by transition between states 
conditional on survival) and the event matrix B. 



Oix(l-T0 
0 
0 



p 

l-p 

1 



0 



0 0 
0 p 
0 0 



l 



n = 



[71 1-71 0] : 



Columns of the matrix IT and <$*¥ correspond respectively to 
state PB, B and dead, columns of the matrix B correspond to the 
observations 'non detected' and 'detected' while rows of the 
matrices fl>*P and B correspond respectively to state PB, B, and 
dead. <I>i denotes the pre-breeder survival, C>2 the adult survival, 
'i'lis the probability of becoming a breeder and p is the detection 
probability. 

Individuals are first captured as PB and then become B; tc 
denotes the proportion of newly marked in state PB. OiX*!*! 
represents the probability that a pre-breeder survives and attempts 
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Table 2. Results of goodness of fit test for the multistate model. 





Test component 


X 2 


d.f. 


P 


Global test 


48.568 


52 


0.610 


Test 3G (test for transience) 


22.264 


23 


0.504 


Test M (test for trap-dependence) 


26.304 


29 


0.609 



doi:1 0.1 371 /journal.pone.0094526.t002 



breeding. For an individual in state B, fl> 2 represents the 
probability that a breeder survives. After the first capture, pre- 
breeders are not observed as they do not attend the colony; 
skippers are not present at the colony; we then estimated breeders 
detection as the detection parameter in the "breeder" state is the 
product of the probability of not skipping breeding and of the 
probability of detection given breeding. 

Goodness-of-fit tests [43] were performed using U-CARE [44]. 
The results of the GOF test are summarised in Table 2. 

We analysed 25 model combinations in which 1* was modeled 
as a constant or functions of age and time, either singly, additively, 
or with interactions. The model assumes that there are two age 
classes (juvenile and adult) and there is an age from which the 
probability to breed ai does not change. Full recruitment age is the 
age at which any individual has become a breeder, it is equivalent 
to age "m" as defined by Gauthier and colleagues [45], which 
considers age "m" as the age at which recruitment probability can 
be considered constant for all subsequent ages. Juvenile survival 
and breeder recapture probabilities were allowed to be constant or 
to vary with time. 

We distinguished the first from the following encounter 
occasions, as the encounter history is conditional on being caught 
in the first period and the following detection probabilities depend 
on the state and the time occasions. We started running a model 
assuming that there is no change after 10 years. As we obtained 
that after age 7 there is no improvement in model fitting we 
reduced the full breeding age to 7 (namely the 6th year of age) 
(Table 3). For a full breeding probability we fixed *P relative to this 
hypothetical age 1. Thus we tested which is the full recruitment 
age (probability of becoming a breeder) by fixing to 1 the potential 
age of full recruitment, specifying initial parameter values we 
constrained beta values before running the model. 

We tested the time effect on survival and the effect of time and 
sex (the last as group effect) on recruitment only. We limited the 
tests in this way because we primarily sexed adults and we have 
information on chicks sex only after 2007. 

Bearing in mind previous results, to assess climatic effects on 
survival and recruitment, model selection began with the model 
with age and time effects on survival and sex effect on recruitment 



age. More in details, after testing the single effect of the individual 
covariates on survival (O) and recruitment (*P) we tested the 
additive effects of different individual covariates and time on 
survival (O) and the additive effects of different individual 
covariates, time and sex on recruitment (*P). 

Models were fitted with the E-SURGE program [46] (for model 
implementation in E-Surge see Material SI). Model selection was 
based on the Akaike Information Criterion (AIC) [47,48]. 

Effects of winter climate on survival over time 

Because data on environmental conditions are available only 
from 2001 to 2010, we used a subset of birds ringed as chicks and 
adults in the same period (2918), modelling potential effect of 
covariates only on adults survival. We carried out goodness-of-fit 
for the single state model using U-CARE [6] . Results of the GOF 
test are summarized in Table 4. Due to the presence of transient 
individuals, we used 2-age-class models and a variance inflation 
factor (c-hat= 1.4). Following the single state, as simplification of 
the scheme reported above, we implemented the model in 
program E-SURGE [44] (for model implementation in E-Surge 
see Material SI). 

We analysed 24 model combinations in which O was modeled 
as a constant or function of time, we then explored singly and 
additively interaction with environmental (external) covariates 
dependent models. The model assumes that there are two age 
classes (young, first year of age, and adult, older than 1 year), 
survival and breeder recapture probabilities were allowed to be 
constant or to vary with time. 

Model selection procedure: model selection was based on the 
Akaike's information criterion (AIC) [47], adjusted for over- 
dispersion (Q). After identifying time-dependent parameters, we 
investigated the effects of SST and CHL in describing the time 
variation with an analysis of deviance (ANODEV) and calculated 
the proportion of variance explained by the environmental 
(external) covariates using the R 2 statistic [35]. 

Ethics Statement 

Field work was carried out under permission from the Marine 
Protected Area. Permits included obtaining DNA samples for 



Table 3. Selection of the full breeding age (f.b.), only best models are presented. 



Model Deviance NP AIC AAIC AIC weight 



f. b. year 7 (age 6) 


3127.74 


29 


3185.74 


0.00 


1.000 


f. b. year 5 (age 4) 


3224.01 


27 


3278.01 


92.27 


0.000 


f. b. year 6 (age 5) 


3224.01 


28 


3280.01 


94.27 


0.000 


f. b. year 8 (age 7) 


3224.01 


30 


3284.01 


98.27 


0.000 


f. b. year 9 (age 8) 


3224.01 


31 


3286.01 


100.27 


0.000 



Deviance (Deviance), number of parameters (NP) and AIC [47] are reported. Preferred model is in bold. 
doi:1 0.1 371 /journal.pone.0094526.t003 
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Table 4. Results of the GOF test for the single-state model. 







Test component 


X 2 


d.f. 


P 


Global test 


233.82 


61 


0 


Test 3.SR (test for transience) 


162.40 


10 


0 


Test 2.CT (test for trap-dependence) 


9.23 


9 


0.42 


Global test (corrected for transience) 


71.42 


51 


0.03 



doi:1 0.1 371 /journal.pone.0094526.t004 



sexing. There were no ethics committee approval requirements to 
carry out this project due to the minimum handling of the 
individuals that was carried out according to Italian and Sicilian 
region legislation (LN 157/92, LR 105/99 and LR 74/2012) by 
certified ringers authorized by the Institute for Research and 
Environmental Protection (ISPRA) to handle adults and chicks. 

Results 

Individual covariates were defined and we made sure that the 
covariates used as integrative climatic indexes were not inter- 
correlated (Table 5). WSR is represented by the first principal 
component explaining 86,49% of the total variance (figure S2a). 
PB is represented by principal components 1 and 2 (PC1PB and 
PC2PB hereafter). The first factor, explaining 42.57% of the total 
variance, is correlated positively with temperature and negatively 
with rainfall, whereas factor 2, explaining 32.19% of the total 
variance, is positively correlated with the NAO (figure S2b). 
During the pre-breeding period, warm and dry conditions are 
associated with a high NAO, whereas cold and wet conditions are 
associated with a low NAO. B is represented by the first principal 
component (PC IB hereafter, explaining 45.97% of the total 
variance), which is positively correlated with rainfall and the NAO. 
B is also represented by the second principal component (PC2B 
hereafter, explaining 33.42% of the total variance; figure S2c), 
which is positively correlated with temperature. These results 
indicate that during the breeding period (i.e., the chick-rearing 
period), warm and dry conditions are associated with a low NAO, 
and cold and wet conditions are associated with a high NAO. 
Trends of individual and environmental covariates are summa- 
rised in figures S3a and S4a. 

By fixing to 1 the potential age of becoming a breeder (full 
recruitment) or the age at which all the birds of the colony may be 
considered breeders, we found that the 2 nd year is the age at which 
most birds start to recruit with a small difference with the 3 rd year 
of age (AIC = 3176. 14 and AIC weight = 0.342; AAIC = 1.60 and 
AIC weight =0.154 for 3 rd year of age, but see Table 6 for more 
details) suggesting that at the 2 nd and, in smaller proportion, at the 
3 rd year of age are Marettimo's storm petrels are becoming 
breeders. In table 6 both the results of full recruitment age 
identification and of sex effect on recruitment age are reported. 
The best model considered constant annual survival trends, 
whereas recruitment was affected by sex, with males recruiting 
earlier than females (see also figure S5). Recruitment probability in 
males is slightly higher than in females during the first 3 years of 
age (in average 4.8%), this difference is no more perceivable from 
the 4 th year. Full breeding age in the Marettimo colony is at the 6 th 
year of age (Table 3). 

We found high mortality in juveniles. The average first year 
survival after fledging was 0.222 (±0.093) and adults survival was 
0.923 (±0.060). Survival was highly variable (figure S3b), 
especially for juveniles. 



Considering single and additive effects of group (sex) and 
individual covariates (local climatic descriptive indexes) we found 
that recruitment differed between sexes with males recruiting 
earlier. Pre-breeding season conditions with colder springs was 
negatively affecting recruitment (Table 7), whereas survival in first 
year after fledging and adults was affected positively by warm 
breeding-period conditions and negatively by abundant rainfall 
explaining the correspondence of the drops in figure S3b with 
higher B and WSR and lower PB values in figure S3a. 

After defining environmental covariates as SST and CHL PCA 
factors (Table 8) we obtained cues on their effects on adults 
survival and encounter probability. We obtained that SST and 
CHL, are not correlated (Spearman R= —0.336, p = 0.321, figure 
S4a) giving strength to the models involving them both. We can 
confirm that SST and CHL may influence storm petrels ecology 
during winter time [39]. A higher sea surface temperature and 
lower chlorophyll concentration were associated with lower 
encounter rates during the subsequent breeding period (figure 
S4b) and a higher encounter rate after favourable years. Winter 
sea conditions in different areas of the Mediterranean affected the 
survival of the Marettimo birds. In particular, we found that 
periods of low productivity (high SST and low CHL) in the 
Alboran Sea negatively influenced adult survival (Table 9). In fact 
we found a significant effect of Alboran Sea conditions 
(AIC = 4175.04 and compared to time effect: AAIC = 1.37, AIC 
weight 0.335, F 2jl5 = 14.385, p = 0.004), with 62% of the temporal 
variance in adult survival explained by the Alboran sea-surface 
temperature and chlorophyll variation. When considering neigh- 
bouring models (SST and CHL covariates considered separately) 
we obtained that SST has a stronger effect (AIC = 41 75.04 and 
compared to time effect: AAIC =1.37, AIC weight 0.335, 
F 2 ,i5= 14.385, p = 0.004, R 2 = 62%) than CHL (AIC = 4183.76 
and compared to time effect: AAIC = 10.09, AIC weight 0.001, 
F 2>15 = 5.560, p = 0.040, R 2 = 38%). Thus SST alone explains a 
large percentage of the total variation on adult survival. This 
finding suggested that the Alboran Sea was the most probable 
wintering site, and that environmental conditions there affect 
overwintering conditions in the Mediterranean storm petrel 
population of Marettimo. 

Discussion 

Recruitment age determination 

Petrels are long-lived organisms with delayed sexual maturity. A 
prolonged pre-recruitment period and extended lifespan provide 
seabirds with sufficient time to explore marine habitats and to 
gather crucial information about prey patches and the spatio- 
temporal variability in their availability [5]. Although we found 
survival values similar to those reported by a previous study [32], 
the recruitment age for the Marettimo colony was lower, perhaps 
due to the greater number of individuals located in this colony and 
the much greater intra-individual variation that resulted. Mar- 



PLOS ONE | www.plosone.org 



5 



April 2014 | Volume 9 | Issue 4 | e94526 



Local and Global Climate Effects on Storm Petrel Demography 







CC 






l/l 












00 


















0. 






03 












(J 






0. 






O 






< 






z 






00 






c 












# c 






*i5 






03 






c 






n 






1- 






00 




Q 






0 






CT1 






.5 

TJ 


PB 




V 


M 




£ 








0. 




0) 






> 






0 






Si 






~* 


00 






0. 






G 






0. 


































o> 














O 




(0 


< 






z 




0 


00 






0. 




3 






























To 


e 




£ 


ro 




o 


c 




u> 


'3 




.2 


cc 




TJ 


CO 




0> 


0. 




.C 






** 






s 






0 


s 




Ol 


a 






1- 






00 






0. 




e 






0 






M 






(5 






Pe 


O 






< 




c 


z 
















a 






c 






0 






ati 






rel 


■3 






J] 




0 


a 




u 


-1 




c 






0 


01 




'■M 


** 




to 


.5 




01 








(5 






> 




0 


0 




u 


u 




on 


tic 




■A 






n 


£ 




HI 






0. 





Ogj 

* I s 
s s a 



o 
< 



fO CL 
<U CD 

2 OU < 
Q. m y 

H Q- £ 



o 
o 
o 



in 
d 



m 

IN 
O 



c! 



U 
0. 







c 

0 






01 _. 




01 _ 




larc 


llatii 


~o 


itic 


din 
pril 


tic 


din 
pril 






sci 


0 


to 

.£ 


$< 


TO 

.£ 




ialy 


~0 

0 


0 

u 


<D 

Q. 




» -5 






E 


<D 








£ S 


01 




0 
c 


Q. 


c 

03 




> 


0. S 


> 


0. S 




ing 


Atl 


II 


rat 


TJ 


rat 




Rainfa 


Breed 
April) 


North 


Pre-Br 
(Marcl 


Integ 


index 
perio 


Integ 


index 
perio 



£ 5 

8 v > 



ra o 

E 'S 



03 
Ol 



o 
o 
o 



^- 
00 
d 



o 
o 
d 



Ol 

o 



o 
d 



5.ES = 

<U 2. ^ 
h CO < 



re en 

= E 

cc CO < 



0 




m 


< 










u 


CO 




0. 


c 






0 








>. 




J2 


re 




'u 




(0 

.£ 


0 

_u 


iod 


"u 




cu 


01 


c 

TO 


Q. 


> 


Atl 


CT1 


rat 


.c 




01 


e 




01 


0 




+j 


z 


£ < 


c 



III 
rM 
U 
0. 



01 01 

.E 3 

■o < 

0> i. 
£ « 
00 £ 



01 Ol 
.£ 3 



03 



£ £ 

_Q ^ IT) 

~ '5 0 

U = O 

1/1 u *^ 

"O °- rs| 

1/1 xs ^1 

0) c J 

is — 

TO o 



01 T3 

S .£ 



™ o 

"5 q 

o E S 

: a 
° -1 

a 0 I 

B .0 

-= u t— 
> S d 

I — U T3 



PLOS ONE I www.plosone.org 



6 



April 2014 I Volume 9 | Issue 4 | e94526 



Local and Global Climate Effects on Storm Petrel Demography 



Table 6. Selection of the full recruitment (f.r.) age and test of time effect on survival and sex effect on recruitment age, only best 
models are presented. 



Model 


Deviance 


NP 


AIC 


AAIC 


AIC weight 




3002.59 


52 


3106.59 


0.00 


0.827 


0, <F 


3011.83 


49 


3109.83 


3.24 


0.164 


4>, *F f.r. year fixed to 3 


3019.74 


48 


3115.74 


9.15 


0.008 


f.r. year 3 (age 2) 


3128.14 


24 


3176.14 


0.00 


0.342 


f.r. year 4 (age 3) 


3127.74 


25 


3177.74 


1.60 


0.154 


f.r. year 5 (age 4) 


3127.74 


26 


3179.74 


3.60 


0.165 


f.r. year 2 (age 1) 


3134.77 


23 


3180.77 


4.63 


0.034 


Deviance (Deviance), nu 


mber of parameters (NP) and AIC [47] are 


reported. Preferred models 


are in bold. 







doi:1 0.1 371 /joumal.pone.0094526.t006 



ettimo colony study conditions and structure are different from 
most of known colonies, being the colony placed in a big and 
completely dark cave accessible only from the sea (thus predator- 
free), and with a large number of accessible nests. The colony 
characteristics make it very different from other monitored 
Mediterranean colonies such as Filfla, Malta, where nests are 
scattered in single inaccessible burrows in the rocks [24], and 
Benidorm, Spain, where the colony is small (ca 300 nests) and in 
two caves easily accessible for predators [49]. The database we 
used in this study is based on a more complete and longer time 
series (2 1 years with one missing year vs 1 5 years with two missing 
years) than that used previously for the same colony in a 
comparative study [32]. This allowed us to investigate and 
describe cues of the subspecies life history, thanks to predators 
absence, making it possible to infer about chicks natural mortality, 
and the long time series allowing correlational analyses with 
climatic conditions trends. The population studied showed that the 
pre-recruitment period may be shortened by unfavourable 
environmental and climatic conditions that may have lead to 
earlier maturation as a consequence of lower adult survival (figure 
S3b). Long-lived birds invest more in survival than fertility [50]. 
However, this study found that storm petrel tend to recruit earlier 
in this colony compared to other monitored colonies throughout 



the Mediterranean [32]. It is possible that this result is due to the 
effects of environmental conditions on survival. 

Effects of environmental conditions on demographic 
traits 

Survival exhibited two troughs (figure S3b), which might be due 
to factors related to colony dynamics (internal) and to environ- 
mental factors (external). The influence of climate on demography 
is represented by the finding that the colony region's climatic 
conditions affect the recruitment age and survival of chicks and 
adults. In particular, chick survival was affected by breeding- 
period conditions and rainfall. Abundant rainfall before and cold 
weather during the chick-rearing period resulted in higher chick 
mortality, as observed in 2010. These conditions may have 
produced wet conditions in the colony caves and reduced 
thermoregulation capabilities, inducing higher metabolic costs 
for chicks and parents [18]. In fact, in 2010, a year with higher 
rainfall and lower temperature, we observed exceptionally high 
chick mortality in the field (only 20 fledged in the surveyed area of 
the colony), while we usually don't find dead chicks when visiting 
the colony. In the view of our results warmer summers and lower 
rainfall predicted as a consequence of the changing climate [51] 



Table 7. Model selection results, only best models are presented. 



Model Deviance NP AIC AAIC AIC weight 



*I*B1+B2+t+WSR *PpB1+PB2+sex 


2984.08 


55 


3094.08 


0.00 


0.842 


®Bl+B2+t ^PBI+PEU+sex 


2989.76 


54 


3097.76 


3.68 


0.134 


®Bl+B2+t ^*PB1+PB2 


2999.35 


51 


3101.35 


7.27 


0.022 




3002.59 


52 


3106.59 


12.51 


0.001 




3011.83 


49 


3109.83 


15.75 


0.000 


®B1+B2 ™ 


3104.72 


29 


3162.72 


68.64 


0.000 


1'PB1+PB2 


3112.46 


30 


3172.46 


78.38 


0.000 




3112.68 


32 


3176.68 


82.60 


0.000 


( I ) PB1+PB2 *P| 


3121.82 


28 


3177.82 


83.74 


0.000 


HWSR f , 


3124.11 


27 


3178.11 


84.03 


0.000 


^WNAO *Pi 


3126.47 


27 


3180.47 


86.39 


0.000 


Hi f , 


3127.74 


29 


3185.74 


91.66 


0.000 



Climatic (PB: Pre-Breeding, B: Breeding, WSR: rainfall, WNAO: winter NAO), "time", and constant effects "i" were considered on survival, <1>, ^skipping the transition 
probability from state pre-breeder to breeder, ¥ ("+" denotes additivity). Deviance (Deviance), number of parameters (NP), AAIC and AIC weight are reported. 
doi:1 0.1 371 /journal.pone.0094526.t007 
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Table 8. Percentage of variance expressed by each principal 
component factor of the PCAs performed (Alboran Sea, 
Balearic Sea, Sicilian Channel, Adriatic Sea, and Aegean Sea: Al, 
L, S, Ad, Ae, respectively; sea surface temperature: SST and 
chlorophyll: CHL). 






°/o COV 


LsstPCI 


70.069 


LchlPCI 


44.937 


LchlPC2 


29.623 


AlsstPCI 


82.34 


AlchlPCI 


86.231 


AdsstPCl 


87.756 


AdchlPCI 


51.557 


AdchlPC2 


29.163 


SsstPCI 


88.377 


SchlPCI 


78.176 


AesstPCI 


88.872 


AechlPCI 


92.015 


doi:1 0.1 371 /journal.pone.0094526.t008 



will not have a negative effect on juvenile survival rate of the 
Mediterranean subspecies. 

Adult survival was directly affected by the climatic conditions at 
the time of breeding due to the adult's needs to meet its own 
energetic requirements as well as those of the chick. Males tend to 
recruit earlier than females and recruitment was reduced in years 
with adverse prebreeding conditions, suggesting that younger birds 
decided to delay recruitment to wait for better condition or more 
productive years. Carryover effects for the winter might occur as a 
consequence [52]. It is noteworthy the fact that males of 
Mediterranean storm petrels are smaller than females [53] and 
consequendy adults may rear the chicks with fewer energy 
expenditure. 

Identification of wintering-area 

The original aspect of our study is the use of statistical models to 
identify seabird wintering areas, as illustrated for the storm petrel, 
a species difficult to follow at sea due to lack of technology, as 

Table 9. Model selection results. 



geolocators or satellite tracking devices are not small enough to be 
used on this species. The trends of sea surface temperature and 
chlorophyll [39] (figure S4a) may influence storm petrel survival. 
In particular, the conditions in the Alboran Sea zone resulted to be 
influencing adult survival. Considering sea surface temperature 
and chlorophyll concentration effects separately, we observed that 
SST has a stronger direct effect while the CHL direct effect is 
weaker (Table 9). It is still unknown what the Mediterranean storm 
petrel feeds on during winter. Our results suggest that they 
probably feed mainly on small fishes, not directly affected by CHL, 
as during the breeding period [27], and not on zooplankton that 
would be directly affected by CHL concentration. Sea conditions 
may induce birds to skip reproduction in years with high surface 
temperatures and low chlorophyll levels at wintering sites, resulting 
in lower encounter rates during the subsequent breeding season. 
This would explain the drops in encounter rate observed in 
breeding seasons after high sea surface temperature and low 
chlorophyll concentration in winter conditions (figure S4a and b). 
Our results highlighted the relationship between environmental 
covariates and encounter rates, as also found in the Pacific and 
other studies on impacts of climate change on avian populations 
[11,54]. 

Finally with respect to predicted features of climate change [51], 
as we observed that colony region conditions affect chick and adult 
survival, whereas conditions in the wintering areas affected adult 
survival, we can conclude that the stability of the Alboran sea 
surface temperature due to anticyclonic systems [55] may result in 
continuous upwelling in the area thus resulting in stability of 
survival of adult storm petrels during winter while the anticipated 
decrease in precipitation and warmer summers [20] will positively 
affect chick survival. Although we did not consider climatic effects 
on the food chain our result suggest future demographic stability 
and an improvement of the species status. Further research on 
climate change effects on the food chain and consequent effects on 
the Mediterranean storm petrel are currently missing and would 
help to forecast the species future conditions. 

Supporting Information 

Figure SI Colony and meteorological stations location. 

(TIF) 

Figure S2 Correlation graphs for principal component 
analysis (PCA) for building integrative climatic indices. 

a) PCA of West Sicily rainfall data; b) PCA of the pre-breeding 





Model 


Deviance 


NP 


AIC 


AAIC 


AIC weight 


time eff 


4125.67 


24 


4173.67 


0.00 


0.664 


Al sst chl 


4145.01 


IS 


4175.04 


1.37 


0.335 


Al sst 


4145.01 


15 


4175.04 


1.37 


0.335 


Al chl 


4156.74 


15 


4183.76 


10.09 


0.001 


L sst chl 


4155.86 


15 


4185.86 


12.19 


0.001 


Ad sst chl 


4161.70 


15 


4191.70 


18.03 


0.000 


Ae sst chl 


4167.40 


15 


4197.40 


23.73 


0.000 


S sst chl 


4170.41 


15 


4200.41 


26.74 


0.000 


Constant 


4175.93 


14 


4203.93 


30.26 


0.000 



Environmental, "time", and constant effects "i" were considered on survival, O, with climatic covariates in 5 areas of the Mediterranean (Al: Alboran, L: Balearic, Ad: 
Adriatic, Ae: Aegean, S: Sicilian Sea). Deviance (Deviance), number of parameters (NP), AAIC and AIC weight are reported. 
doi:1 0.1 371 /journal.pone.0094526.t009 
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climatic variables (temperature, rainfall, and NAO); c) PCA of the 
breeding season climatic variables (temperature, rainfall, and 
NAO). Abbreviations are as in Table 6. 
(TIF) 

Figure S3 Trends of individually based covariates and 
demographic traits, a) Time series of climatic indices 
considered as candidate covariates in the analysis of survival and 
recruitment (PB: Pre-Breeding, B: Breeding, WSR: rainfall, 
WNAO: winter NAO). b). Estimated survival at Marettimo, 
confidence intervals in dotted lines. 
(TIF) 

Figure S4 Trends of environmental covariates and 
demographic traits, a) Time series of sea surface temperature, 
SST, and chlorophyll concentration, CHL, in the Alboran Sea. b) 
Estimated survival and encounter rate at Marettimo's colony using 
Alboran Sea climatic conditions (CI values are nil thus not visible 
in the graph). 
(TIF) 
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